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Abstract 

A semianalytic method to estimate the angular resolution of tracks, that have been 
f*-. I reconstructed by a likelihood approach, is presented. The optimal choice of coordi- 

. nate systems and resolution parameters, as well as tests of the method are discussed 

^ . based on an application for a neutrino telescope. 
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1 Introduction 



This paper describes a statistical procedure to extract resolution estimates on 
a track by track basis. The method was developed for the AMANDA Neutrino 
Telescope at the South Pole [3], which uses a 3-dimensional grid of photosen- 
sors imbedded in highly transparent ice to provide spatial and time resolution 
of Cherenkov photons, that e.g. arise from long muon tracks. 

The knowledge of track resolutions is of particular importance in the search for 
localized sources, such as distant galaxies. The resolution information can in 
addition be used to suppress mis-reconstructed tracks, that typically are less 
well defined. The method is not limited to muon reconstruction in neutrino 
telescopes and can be applied to any experiment in which tracks have been 
reconstructed with a likelihood approach. 
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2 Technical aspects of obtaining the confidence ellipse 



2. 1 Definitions and prerequisites 



In all following paragraphs a pattern of hits in an event will be modelled by an 
infinitely long oriented straight line. It will be parameterized by an arbitrary 
point r on the track and the zenith i? e [0, tt] and azimuth 9? e [0, 27r] angles 
in a given detector coordinate system. 

The procedure presented operates on likehhood functions, thus a hkelihood 
based reconstruction must be available. For any given hypothesis {ro, t^o, <^o} 
and an observed pattern of hits P, a likelihood £(ro, "(^O) V'o; P) must be cal- 
culable. The hypothesis hbest = {i"best, ''^best, V'best} with the highest value of 
C = Chest is considered to be the best estimate for the true values. It is as- 
sumed that this extremum has already been found. 

The uncertainties, which are to be determined, are in essence ^ the errors of 
■J? and (fi. It is mandatory to use the optimal coordinate system for each event 
and to consider the errors in the determination of r as well. 

In order to arrive at a confidence ellipse in the two dimensions of track direc- 
tion, the set in parameter space is being searched for, where the value of the 
negative logarithmic likelihood function has changed by 1/2 with respect to 
-log£best: 



A(-log£) = (-log£ellipse) - (-logAest) = 2 • (1) 



After discussing the coordinate systems used, the dimension of the problem 
will be reduced to the two dimensions of direction by eliminating the three 
dimensions of the location of r. This is done using a numerical minimization 
process. The remaining ^^-minimization is done analytically to obtain the 
uncertainty estimators. 

The shape of the negative log-likelihood function around its minimum is con- 
sidered Gaussian. The tests described in section 4 prove that this assumption 
is realistic. 
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Figure 1. Figure a) shows a track in a Cartesian coordinate system {x,y,z). The 
direction towards the track's origin is described with standard spherical coordinates 
("i?, (f). The downward orientation of the z-axis follows the AMANDA standards 0. 
In b) a rotated coordinate system (x', y', z') is presented, where the x'-axis is defined 
by the track. To describe directions that are close to the track's direction, relative 
coordinates 9 = — 11/2 and (p = ip' — ti are introduced. They are approximately 
the coordinates in the tangent plane, as shown in c). 

2.2 On coordinate systems 

There are two sets of coordinate systems used in the method presented here. 
The first system uses all five coordinates {x,y, z,{},ip), which are connected 
to the 3-dimensional Cartesian detector coordinate system as illustrated in 
Figure 1 a). In the second system the problem is reduced to the 2-dimensional 
subset of angles and 6. 

The complication with i} and if is that the first has boundaries and the second 
is periodic. Moreover, at the boundaries of the variable if becomes meaning- 
less. In a mathematical sense, i) and cp are locally orthogonal. However, if one 



A coordinate transformation needs to be done first, which will be explained in 
the following section. 
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looks at these coordinates in the area surrounding a point, they are distorted 
- in particular close to -d — and -d — tt. 

As the behaviour of the likelihood function around the well known best track 
estimate hbest is to be investigated, the coordinates -d and ip should be as close 
to being orthogonal as possible at that particular direction ("i^bcstj V^bcst)- Hence 
the 3-dimensional detector coordinate system is rotated such that the new x- 
axis, x', is defined by the track. This results in t^^gg^ = | and v^^est — ^ 
illustrated in Figure 1 b), where the distortion is minimal. 

Here and in the following section 3 only relative coordinates 9 = d' — tt and 
(j) = ip' — tt/2 with respect to t^^gg,. and <^best ^^^^ ^e referred to. This can be 
thought of as local coordinates in the tangential plane to the unit sphere at 
i'&best, fhest) with the Origin at the point, where the plane touches the sphere ^ . 
By choosing a proper rotation, it is assured that the ^-direction is parallel to 
the ■(9-coordinate, i.e. towards the negative z/z'-a.xis. The rotation must be 
such that the track, the 2;-axis, and the ^'-axis are in the same plane. 

The local coordinates are displayed in Figure 1 c). 

2.3 Reducing the dimensionality of the problem 

The location of the point r is of no concern for the discussion of where a 
track points to. Nevertheless the correlations of the errors in direction with 
the errors in determining r need to be considered. 

The reduction procedure of the dimension shall first be discussed assuming 
only one spatial dimension x and one directional dimension 9. 

Consider the two-dimensional confidence ellipse in (^,a;)-space (Figure 2). It 
is tilted versus the ^-axis, which indicates a correlation of 9 and x. For each 
9 one determines the value of x for which the common likelihood function 
gets maximal. In the strictly Gaussian case, the result is a straight line. The 
likelihood function along this line is again a Gaussian as a function of 9. The 
width of this Gaussian is ae as can be seen by straight forward calculation. In 
this error estimator ae the correlation cov{9,x) is included by construction. 

The extension to three spatial and two directional dimensions is straight for- 
ward: first one determines a set of points in (0j,6'i). For each point, the best 
guess, i.e. the best hkelihood k = — log JCi, is found with respect to r. 

^ This is an approximation that imphes that the coordinates on the surface of a 
sphere can be assumed Cartesian for small distances. The maximum error intro- 
duced is of the order of 0.15% for distances less than 5°. 
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Figure 2. Confidence ellipse in two dimensions with a directional coordinate 9 and 
a spatial coordinate x. The rectangle around the ellipse has the size 2 erg by 2ax- 
Finding for each the corresponding x with extremal likelihood, one arrives at the 
line shown. 

In the Gaussian approximation these points li{0i.9i) form a paraboloid. Its 
parameters are found by analytic minimization, which is discussed next. 



2.4 The analytic minimization 



In order to find the paraboloid that best fits the points li{(f)i, 6i), an analytic x^- 
minimization is performed. The paraboloid 9; 6, P2, T) is parameterized 
in its polynomial representation: 

6; 5, A, r)^5 + M + M + ^(0, ^)r(0, ef . (2) 
The constant parameter is represented by S and the parameters in first order 



contains the second 



by Pi and ^2- The symmetric 2 x 2- matrix F = 

^ri2 T22 f 

order parameters. Using the method of least squares one calculates 

n n 1 

1=1 i=i ^ 

(3) 

and requires the minimization conditions 

^ ^ ^ = ^ L n iw 
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This leads to a six-dimensional inhomogeneous system of linear ^ equations: 



E0i^i iE0f E'/'f^i I EM' 
E ^1 iE0?^^ EM' |E 

E<P^O^ E<P10, E(|>^0^ ^E^PP, E(l>l0f lE(|>^Of 
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(5) 



One can choose the set of points (0,, 6i) in a suitable way to simplify equa- 
tion (5). If for each point {4>i,0i) also {(f)i,—9i), {—(f)i,9i), and {—(f)i,—9i) are 
added to the set ^ , one arrives at 
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(6) 



This leads to three equations: 

j29lP2=T.h0i 

J2</^^9l{^U2 = T.h<P^0^ 

and a three-dimensional system 



(7) 
(8) 

(9) 



^ These equations are linear in the parameters S, P2, and F, which are to be 
determined. 

^ This is both a sensible and easily fulfilled demand, as h^^^^ corresponds to ^ = 
(t) = 0. 
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Figure 3. The confidence ellipse in 6 and (j) can be represented either by ag, a^, and 
the covariance cov{9, (p) or by the major and minor axes ai and (T2, and the rotation 
angle a of the ellipse between the major axis and the ^-direction. 
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which can easily be solved. 

The error estimators sought are contained in the parameters of the paraboloid. 
The three degrees of freedom of 6 and /?i,2 correspond to the position of the 
minimum of the curve. The three parameters in F contain the curvature of the 
paraboloid, which are functions of the resolution, as the covariance matrix C 
of the Gaussian approximation is the inverse of F. From C = F~^ one obtains: 



cov 



(r-^)ii 
(r-^)22 
(r-^)i2 



(11) 

(12) 
(13) 



The results can also be represented in a more intuitive way using the major and 
the minor axes ai and cr2 of the confidence ellipse, and its rotation angle a with 
respect to the 0-axis (Figure 3). In order to find this alternative representation, 
the covariance matrix C needs to be diagonalized. That leads to 



cr 



1,2 




+ detC 



(14) 
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Choosing ai > (72, the rotation angle a is then given by the expression 



tana = 




(15) 



In the following discussion a third representation (ua, e, a) will be used, where 
(Ta is connected to the area of the ellipse and e is its excentricity: 



The rotation angle a is the same as in the previous parameterization. 

This procedure can lead to unphysical, negative values for the variance es- 
timators. This happens whenever the basic track finding procedure has not 
arrived at a good minimum in likelihood space. The search for the minimum 
should then be repeated more carefully. 

Note, that the resolution estimates are calculated for each pattern of hits 
seperately. Still they are not the true values for the discrepancies of true and 
reconstructed direction, but estimates of the average discrepancy for a set of 
events, which are similar to the one under consideration. 



3 Space angle resolution 

For a neutrino telescope it is common to quote the median ^ of the spatial 
angle ^ between true and reconstructed track direction in order to describe 

its resolution. This section is devoted to the connection of the error ellipse 
parameters to the space angle uncertainty ip. It is convenient to choose the 
representation of the ellipse using Ua and e, with a not playing a role in this 
special context. 

The true direction lies within the area given by the error ellipse with a proba- 
bility of 39.4%. The median of a distribution corresponds to a 50%-probability. 
Hence as a first step the ellipses need to be enlarged correspondingly. 

In order to obtain ip, one considers the Gaussian likelihood function in its 




(16) 
(17) 
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Figure 4. Obtaining the space angle resolution estimate by applying a correction 
based on the excentricity. For fixed (Ja = 1° and varying e, the equivalent radius 
M is displayed, such that the true direction lies within a circle of radius M with a 
probability of 50%. 

diagonal form, where the coordinates^ are xi and X2: 

To obtain the median M, the function f{xi,X2) must be integrated from the 
center (0, 0) outwards until a cumulated probability of 0.5 is reached. This 
integral is done in polar coordinates (radius p and angle 7) 

0.5 = / / / (pcos7, psin7) pd7dp . (19) 
Jo Jo 

For the special case of o"a = (Ti = (T2 (i.e. e = 1), the calculation can be carried 
out analytically and yields M = 1.177 • Ua. 

For the general case (with arbitrary e) the integration must be carried out 
numerically. This has been done for the case of Ua = 1° and the result is 
shown in Figure 4. The median M scales linearly with (Ja. 

The curve is symmetric with respect to the exchange of (T\ and (72- To im- 
prove the speed of the program an approximation proportional to (e + ^) has 
been applied which results in the excentricity corrected spatial angle resolution 
estimation: 

< = 0.57-(e + V^a ■ (20) 

The fit is also shown in Figure 4. The deviation in the interval e G [1;3] is 
below 5%. For applications where this is too coarse an approximation, the 
exact results of the numerical integration can be stored in appropriate tables. 

^ These coordinates xi and X2 are connected to cj) and 6* by a rotation of the angle 
a. 
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4 Discussion 



The above algorithms have been tested within the framework of the AMANDA 
neutrino telescope. The estimates perform in a stable and reliable way and 
produce sensible results. 

To ensure the correctness of the estimations, several testing procedures both 
in data and in Monte Carlo have been carried out: 

(1) For (Tg and a^f, its pull has been studied. It is defined as the ratio of the 
difference of true and reconstructed direction over the resolution estima- 
tor. If the estimation is good, the pull should be Gaussian distributed, 
centered at zero and with unit width. This investigation can only be done 
in a Monte Carlo simulation. 

(2) For the spatial angle resolution the pull is not a sensible quantity, because 
angles in space can only be positive. Hence one can study ensembles 
of events with the same cXa or a^, respectively. The median of the true 
deviation can then be compared to the corresponding a^,. Again this can 
only be done in the simulation. 

(3) In data the true directions are not known. Thus another approach is being 
followed. Each event is split into two subevents by assigning every second 
hit to subevent 1 and the remaining hits to subevent 2. Each subevent 
undergoes reconstruction and is subjected to the resolution estimation 
algorithm. That leads to directions^ t?*, y?* and errors cr^ und cr^ with 
ZG{1,2}. 

The pull 

P„ = ^ = , 21 



should also be a Gaussian distribution centered at zero with unit width. 
For the difference in azimuth angles must be multiplied by sin'd to 
make up for the differences in the rotated and not rotated coordinate 
systems: 



Of course this last check can also be done in the simulation. 

In all tests the results were satisfactory. Further details about the tests and 
their results can be found in 



Note that the reconstructed directions are expressed in standard detector coordi- 
nates, whereas the error estimates are obtained in their respective rotated systems. 
That is due to technical reasons only. 
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4-1 Concluding remarks 

(1) Often one applies quality cuts to the data set, which rely on additional 
information that is not reflected in the construction of the likelihood 
function itself. Hence this information is also not used in obtaining the 
resolution parameters. In this case one expects the selected tracks to be 
on average closer to the true direction than from the likelihood analysis 
alone. Consequently one observes a pull with a width smaller than 1. 

(2) The resolution parameters - such as Ua - can efficiently be used as cut 
variables, as misreconstructed events on average have a worse resolution 
than correctly reconstructed ones. 

(3) The resolution estimation on an event-per-event basis can be used as 
additional input in point source searches. The information is hard to be 
included in the commonly used binned search algorithms [2j . An appropri- 
ate procedure based on maximum likelihood methods that can integrate 
the new information in a natural way will be presented in a forthcoming 
publication. 
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